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We analytically study the phase transitions between s-wave holographic insulator/superconductor 
and metal/superconductor. The problem is solved by the variational method for the Sturm-Liouville 
eigenvalue problem in the theory with dark matter sector of U(l)-gauge field coupled to the Maxwell 
field. Additionally in the probe limit we investigate the marginally stable modes of scalar pertur¬ 
bations in the AdS solitonic background, connected with magnetic field in the dark matter sector. 

We have found that even with dark matter sector the superconducting transition temperature Tc is 
proportional to charge density p in power 1/3. This value seem to be strong coupling modification 
of the exponent 2/3 known from the Bose - Einstein condensation of charged local pair bosons 
in narrow band superconductors. The holographic droplet solution is affected by the coupling to 
the dark matter. Interestingly in the probe limit the critical chemical potential increases with the 
decreasing coupling to dark matter making the condensation transition harder to appear. 

PACS numbers: 11.25.Tq, 04.50.-h, 98.80.Cq 


I. INTRODUCTION 

The gauge/gravity duality in the form of AdS/CFT correspondence [Il,i provides an interesting framework to 
study strong coupling effects in quantum many body d-dimensional systems [3| by means of the d -I- 1 dimensional 
spacetime with the negative cosmological constant. In particular, this technique has been widely used to describe phase 
transitions from the normal or insulating to superconducting state. In the original work on ’building the holographic 
superconductor’ ii, the description of the single band s-wave superconductor was proposed. The scalar complex 
field with the appropriate potential has been incorporated into the theory of gravity and the condensation of its dual 
operator at finite temperature T, lower than the critical one Tc was observed. The condensed operator has been 
identified with the superconducting order parameter. In the bulk, the temperature was introduced as the Hawking’s 
black hole one. 

The aforementioned approach has been extended in many directions, taking into account the relevant aspects of the 
existing superconducting materials. For instance, the models of d-wave [1] superconductivity have been considered, 
to shed some light on the strong coupling behaviour of the well known high temperature superconductors J7| , which 
feature this symmetry of the order parameter. The spin triplet superconducting states of simple p-wave [8| as well 
as the chiral Px ± ipy symmetry have also been elaborated in considerable details Q. In Refs.p^ the multi-band 
superconducting systems have been studied in view of many materials in which the coexistence of different orbitals 
plays a crucial role, e.g., in MgB 2 [HI, Sr 2 Ru 04 [l^ or in heavy fermion superconductors [l^. The description of 
these materials requires at least two hybridized orbitals which in the gravity approach translate into two scalar fields. 

On the other hand, apart from the aforementioned studies of conductor/superconductor phase transitions, the 
holographic insulator/superconductor transitions attracted a great attention. Modification of the bulk gravity theory 
by considering the five-dimensional AdS soliton line element [l^ coupled to Maxwell gauge field and scalar one, 
allows for building a model of holographic insulator/superconductor phase transition at zero temperature [T^ . In 
gauge/gravity duality description, the AdS soliton is dual to a confined field theory with a mass gap, mimicking 
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an insulator phase [l^ . It was revealed that in the presence of a chemical potential in the solitonic background, 
the insulator/metal transition is of the second order. Namely, for the chemical potential greater than some critical 
value, the considered background turned out to be unstable and non-trivial hair emerged. This fact is interpreted as 
insulator/superconductor phase transitions. 

In Refs. [Oj 13 it was shown that the strength of various kinds of matter backreactions could generate new types of 
phase transitions. Marginally stable modes of scalar/vector perturbations in the AdS solitonic spacetime were studied 
in [3j 113 to reveal the onset of the phase transition as well as to hnd the magnetic field effects on them. Among 
others, it was claimed that magnetic held made the phase transition harder to occur. The compatibility with the 
earlier investigations were also announced. 

Recently, the inhuence of nonlinear electrodynamics on the holographic insulator/superconductor phase transitions 
was taken into account (^ . , while the problem of p-wave symmetry of the transitions in question was treated in 

Ref. [ 2 ^, [23 • On the other hand, the analytical investigations tackling the phase transitions of this type, in Gauss- 
Bonnet gravity were discussed in 

Moreover, superconducting solutions in which the condensate was conhned to a hnite region and decayed rapidly 
outside during conductor/ superconductor phase transitions, were examined. Both the vortex and the droplet models 
were constructed for s-wave type of superconductors A non-Abelian droplet solutions emerging during insula¬ 
tor/superconductor phase transition in p-wave and p-|-ip-wave symmetry were studied in Ref. [13 . It has been shown 
that in the case of Gauss-Bonnet background the coupling constant of the theory (like in 113) affects the transition 
in question. 

Furthermore, the question of the possible matter configurations naturally appears in AdS spacetime. The problem 
of the strictly stationary Einstein-Maxwell spacetime with negative cosmological constant was treated in [SOjj while 
the simply connected Einstein-Maxwell-axion-dilaton spacetime with negative cosmological constant and arbitrary 
number of C/(l)-gauge fields was examined in Ref.[3l|- It was revealed that the considered spacetimes could not allow 
for the existence of nontrivial configurations of complex scalar fields or form helds. 

Motivated by the above problems as well as to provide continuity with our previous studies [H, [13 , we address here 
the problem of phase transitions among insulator/superconductor and metal/superconductor for s-wave holographic 
superconductors in the theory in which dark matter sector is coupled to the Standard Model. We shall look for the 
imprints of dark matter sector in possible holographic experiments. 

The importance of examinations of such kind of models goes back to the need of explanation of 511 keV gamma 
rays astrophysical observations made by Integral/SPI [13 as well as the experiments showing the electron positron 
excess in galaxy, revealed by ATIG/PAMELA [35|, 113 • Their energies vary from a few GeV to a few TeV depending 
on the experiments. On the other hand, the new physics can explain the 3.her discrepancy between measured value of 
the muon anomalous magnetic moment and its prediction in the Standard Model [37|- The other facet concerns the 
fact that dark matter model is subject to the key ingredient in early Universe, where the topological phase transition, 
giving rise to various topological defects, might have happened. 

Our analysis will be addressed to the theory in which, apart from the gravitational action given by 

^9 = 1 ^ 

where = SttGs is five-dimensional gravitational constant, A = —tSjl? stands for the cosmological constant, while 
L is the radius of the AdS spacetime, we shall examine the Abelian-Higgs sector coupled to the second C/(l)-gauge 
field 

Sm = J V^d^x(^- - vw - , ( 2 ) 

where the scalar field potential satisfies U(t/) = -I- = 2V[^Ai,] stands for the ordinary Maxwell 

field strength tensor, while the second C/(l)-gauge field B^^, is given by = 2V[^B^y Moreover, m, q represent a 
mass and a charge related to the scalar field ip. Here a is a coupling constant between C/(l) fields. The compatibility 
with the current observations establishes its order to 10“^. 

Within the above model the backreaction problems of the dark matter sector on s-wave holographic superconductor 
was analyzed in Ref. [13 • It was revealed that the dark matter coupling constant is bigger the smaller is the critical 
temperature. The so-called retrograde condensation takes place for the negative value of the aforementioned constant. 
In Ref. [13 the nature of the condensate in external magnetic field and the behaviour of the critical field near the 
transition temperature were examined. The obtained upturn of the critical field constitutes the fingerprint of the 
strong coupling. In that study a has been found to be limited to positive values. 

The organization of the paper is as follows. In Sec.II we start by studying the s-wave holographic zero temperature 
insulator/superconductor phase transition using the solitonic AdS background. The chemical potential is the control 
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parameter of this phase transition. In Sec.Ill, black hole background taken as the gravity configuration allows for the 
analysis of the transition observed for // > /i^ from holographic metal at high temperatures (T > T^) to holographic 
superconductor at temperatures below T^. The effect of the dark matter sector on the insulator/superconductor 
transition of the droplet is studied in Sec.IV. We end the paper with summary and discussion of the obtained results 
in the light of variety of the existing metal/insulator and metal/superconductor transitions in condensed matter 
systems. We limit our analysis of the phase transitions to the probe limit. 


II. PHASE TRANSITION INSULATOR/SUPERCONDUCTOR 

In this section we analyze the model of s-wave holographic insulator/superconductor phase transition in five¬ 
dimensional spacetime, where the matter sector is coupled to another C/(l)-gauge field, representing the dark matter 
sector. In the probe limit, we setup the considered model in the AdS soliton background 0 , which line element is 
subject to the relation 


f{r) 

where /(r) = — r^/r^. The geometry resembles a cigar, if one gets rid of (r, (/?)-coordinates, with a tip located at 

r = tq. The AdS soliton solution is achieved by making two Wick rotations on a five-dimensional AdS Schwarzschild 
black hole line element. The asymptotic AdS spacetime tends to x topology near the boundary. A conical 
singularity at cq can be removed by the Scherk-Schwarz transformation of (/^-coordinate, i.e., (/j ~ (/j -|- tt L/rg. Due 
to the compactification of (^-direction, the AdS solitonic background alllows for a description of a three-dimensional 
field theory with a mass gap, which echoes an insulator in the condensed matter physics. The temperature in the 
solitonic background is equal to zero. 

Without loss of generality we put L = 1 and for simplicity we assume that At = 4>(r), Bt = r]{r) and ip = 'ip{r). 
The underlying system of differential equations for scalar and gauge fields yields 
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( 5 ) 
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~ 2 

where we set 5=1 — ^ and ci as an integration constant. 

Next we impose the boundary conditions on the adequate quantities. Namely, at the tip of the AdS soliton we 
demand that the solutions will be provided by 


0 = ipo +-piir - ro) + ■p2{r - rof + ..., (7) 

(p = (po + (pi{r - ro) + p 2 {'r - rof + ■ ■ ■, ( 8 ) 


where 'pm and pm, for the range m = 0, 1, 2,..., are constants. Moreover, in order to achieve the finiteness of the 
considered quantities, one has to fulfill the Neumann-like boundary conditions (0i = 0 and pi =0). Contrary to the 
AdS-black hole case, where at the event horizon p is equal to zero, here it can acquire a non-zero value at the tip of 
the AdS soliton. On the other hand, near r —>■ oo, we have the following behaviours: 


'P_ 

X- ^A+ ’ 


</> = M 


P_ 
2 ’ 


( 9 ) 


where /i and p stand for the chemical potential and charge density in the dual theory, while A± = 2 ± ^4-1- m^. The 
coefficients are responsible for the vacuum expectation values of the operators < 0± > dual to the scalar field. 
One can impose the conditions that either p~ or 'p~^ vanishes [s^. In what follows we shall assume that p~ vanishes 
and consider 0+ =< Oi > with < Oi > denoting the expectation value of the corresponding CFT operator. 

It will be convenient to rewrite the above equations in terms oi z = rg/r variable. They reduce to the forms 
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where the prime denotes the derivation with respect to z-coordinate. 


( 11 ) 

( 12 ) 


A. Critical chemical potential 

It was revealed in Ref.[l^ that when the chemical potential exceeds a critical value, the condensation will set in. 
This state can be interpreted as a superconductor phase. In the case when /i < the scalar field '0 achieves value 
close to zero and the phase can be interpreted as the insulator. The system has a mass gap, which is connected with 
the confinement in (2 + l)-dimensional gauge theory via performing the Scherk-Schwarz compactification. In the light 
of these facts, the critical value of the chemical potential is the turning point in a superconductor phase transition. 

For the chemical potential /x = /Xc, the scalar field is very small ip ^ 0 and the equation (EH) for the gauge field (p 
reduces to the form 

<(>" + (7 + 

The general solution of this equation can be easily found to read 

(p{z) = di + dz log I , ^2 , (14) 

1 + 


where di and ^2 are integration constants. 

In order to fulfill the assumed boundary conditions ([8]) at the tip z = 1, we require ^2 = 0. Thus (p has the constant 
value /X, when ip{z) = 0. Moreover, from the equation ([5]), one obtains that in the considered case p = 0. These 
results are in accord with the numerical analysis presented in flB |. 

By virtue of the above, as /x —>■ /Xc, we have 


r 






f z^ f 


Ip = 0. 


(15) 


By introducing a trial function near the boundary z = 0, in the form ip{z) =< Oi > z^'F{z), where x = + or —, 
and by imposing the boundary conditions F(0) = 1 and F'{0) = 0, the underlying equation can be brought to the 
following form 


(p(z) F'(z))' - q{z) F{z) + 7 r(z) F{z) = 0, 


(16) 


where the various terms in the above relation are provided by 


p{z) 

q{z) 

r{z) 


_^2Ai-2 f ^ii^i ~ 1 ) / 
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(17) 

(18) 
(19) 


According to the Sturm-Liouville eigenvalue problem, we can specify /x^ as a spectral parameter and estimate its 
minimum eigenvalue by varying the following functional 


, J,^dz[F'izrpiz) + qiz)F\z)] 

f^^dzr(z)F^z) 


( 20 ) 


The trial function will be set in the form F(z) = 1 — az^. Importantly, the critical value of the chemical potential 
is unaffected by the dark matter sector parameters. The value of the /x^ in equation (1201) depends on the parameter 
a entering the trial function. Changing a we find numerically minimal value of for a = amin- Both Omm and 
pPp = p^{amin) depend on the parameters of the model - in particular m^. In Table U the critical chemical potential 
has been presented for a few typical values of fulfilling the Breitenlohner-Freedman bound = 
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— 16/4 required for the stability of the AdS^+i spacetime. The same result was analytically obtained in s-wave 
holographic/superconductor phase transition studies in Einstein-Maxwell scalar theory [^ . and is in accord with 
the numerical examinations provided in Ref. M- On the other hand, in Gauss-Bonnet gravity, for the transition in 
question, one observes the critical potential increase with the growth of curvature corrections, for the same mass of 
scalar field. For the fixed value of the strength of curvature corrections, with the increase of scalar field mass, the 
critical potential becomes larger 


A+ 



O^min 

5 

2 

15 

4 

1.890 

0.330 

6 

2 

12 

4 

2.398 

0.371 

7 

2 

7 

4 

2.903 

0.407 

8 

_2_ 

0 

3.406 

0.439 


TABLE I: Values of the critical chemical potential together with a parameters minimizing the functional (1201) . In all 
the above examples we put g = 1.0 and ro = 1. The numerical values of m? and A+ are chosen for illustration 
purposes. They obey the physical requirement that the masses do fulfill the Breitenlohner-Freedman bound 
> —d^/4 = —16/4 required for the stability of the AdS^+i spacetime. 


B. Critical phenomena 


In this subsection, we shall concentrate on studies of the critical exponent for condensation operator as well as on 
the mutual relations between the charge density p and the chemical potential. The question we are asking here is how 
the order parameter of the superconductor i.e. < Oi > and the charge density p depend on the distance [p, — p^)- 
Having in mind the form of the scalar field (V'(-z) =< Oi > z^'F{z), where F{z) is the trial function introduced 
earlier) near the boundary 2 = 0, when p —>■ pc, the relation for gauge At = ^(r) field can be rewritten as 

</" + f 4 + -) '^ = 0- (21) 

\J zj af 

To proceed further, let us remind that for p slightly above the critical value the condensation scalar operator < Oi > 
is very small. This enables us to seek the solution in the form 

<p{z) ^ Pc + <Oi> x{z) + ■■■ (22) 


In order to recover the previous result 4>{z) = p, we have to impose boundary condition x(l) = 0. On the other hand, 
close to the boundary z = 0 one expands the function x{z) = x(0) -I- x'{0)z + ^x"i^)z^ + ■ ■ ■, rewrites the relation (|9]) 
as 


(j){z) ~ p- p z^ cz pc+ <0^> ^x(O) -k x'(0)z -k ^x"(0) + • •. ^ , 

Comparing the coefficients of the z^ and z^-terms, in the above equation, one obtains the relations 

p- Pc~ <0,> x(0) 

x'(o) = o. 


Inserting the relation (|2^ into (I2IL one can easily find that x(z) will satisfy the following equation: 

X"(z) + f 4 + -) ^'(^) - ^(^) + ■ • ■ = 0- 

\fzj a f a f 


(23) 


(24) 

(25) 


(26) 


Close to Pc the term quadratic in < > is much smaller than the linear one and may be safely neglected leading to 

x"(^)+f4+-) ^'(^)=^ ■ 

V/ zj af 


( 27 ) 
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In the next step, let us redefine xi^) function by the new one ^{z) multiplied by the adequate factor 

x{z) = 2 e(^)- ( 28 ) 

a 

It remains to be checked if the new definition of x{z) enables us to get rid of the a in the equation (|27)) and to extract 
the quantity < Oi >. It can be inspected that ^(z) will satisfy the following relation 

+ F^z) = 0, (29) 

and consequently the scalar operators < Oi > imply 


<o, >= 


'{fi - He) d 


2 He m ■ 

In order to find ^(0) it will be helpful to rewrite the equation (12^ in the form 


fzej = 

Then, having in mind the fact that C^(l) =0, leads to the conclusion that ^(0) is provided by 

C(0) = Cl - ^ (^C 2 + dyq^ rl ^^(y)^ . 


(30) 


( 31 ) 


(32) 


The integration constant Ci, C 2 are determined by the boundary conditions imposed on y(z)-function. We relegate 
their determination to the Appendix A. 

On the other hand, the above relations reveal that the operators < Oi > yield 

< C>j >~ T{h-Hc)K (33) 


where the T factor contains information of the dependence on dark matter sector. The bigger a (the smaller value of 
the a-coupling constant we take) we consider, the greater factor one obtains. 

Moreover, our analytical results show that the holographic s-wave insulator/superconductor phase transition rep¬ 
resents the second order phase transition, with the critical exponent of the system attaining the mean-field value 
1/2. The same conclusions were achieved in the case of of the ordinary s-wave holographic insulator/superconductor 
phase transition studies in Refs. Eng. On the other hand, the same form of the dependence was also obtained in 
Gauss-Bonnet theory [23| , confirming the previous numerical results [1^ . The Gauss-Bonnet coupling constant 

connected with the influence of curvature corrections, does enter in the multiplier factor of the scalar operators, but 
the critical value of the exponent takes the mean-field value. 

Next, we find the dependence of the charge density p on the critical chemical potential. In order to calculate p, 
we use (PSII which implies ^'(0) = 0 together with the previous requirement ^(1) = 0 being subject to the boundary 
condition. Comparison of the adequate coefficients of z^-order in equation (1^ . gives the relation for the charge 
density in the following form 

P=-^^X"(0). (34) 

In order to find x^^(O) we rewrite equation (1271) in the form which implies 

(/^X') - F\z)He<0,>=Q. (35) 


Integrating both sides of it and taking into account the aforementioned boundary conditions, one obtains 

^"(0) = ^ 1,^0=-2 < O, > ^ [ dz F\z). (36) 

Then, by virtue of the equations (1361) and (I34L having in mind the relation (I30|) . one arrives at 


P= (m-Mc) b, 


( 37 ) 
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where the quantity B yields 


B = 


2 5 ( 0 ) 


dz z 


2Ai-3 


F\z). 


(38) 


In view of the above relations, the charge density is proportional to the difference Fi (/r — fj,c), where the factor is 
independent on d characterizing dark matter sector and accomplishes the ordinary dependence of the form p ^ {p— 
achieved analytically in [2^ and numerically in Refs. [13) Some typical values of B factor are presented in Table 
IE The dependence of the T factor on the dark matter sector coupling a is depicted in FiglTJ 

In the case of Gauss-Bonnet gravity sector, the factor standing in front of (/i — md, is function of mass and 
Gauss-Bonnet coupling constant. However, the form of the linear dependence survives [39l|. 


A+ 


Me 

5(0) 

r 

B 

5 

2 

15 

4 

1.890 

0.081 

1.801 

1.329 

6 

2 

12 

4 

2.398 

0.062 

1.823 

1.144 

7 

2 

7 

4 

2.903 

0.049 

1.863 

1.029 

8 

2 

0 

3.406 

0.040 

1.913 

0.948 


TABLE II: Values of the prefactors for the condensate and charge density. In all above examples q = 1.0 and trial 

function was of the form E = 1 — az'^. 


III. HOLOGRAPHIC METAL/SUPERCONDUCTOR PHASE TRANSITION 


In this section we shall scrutinize the problem of holographic s-wave metal/superconductor phase transition at low 
temperatures provided by the black hole background. The problem of s-wave holographic superconductor with dark 
matter sector in the context of the backreaction of matter fields on the gravitational background was investigated in 
[33 |. where the critical temperature was found. On the other hand, for n-dimensional gravitational background, the 
expectational value of the scalar operator and the influence of magnetic field on the holographic superconductor was 
analyzed [s^. For the completeness of the investigations, using quite different methods, we elaborate the dependence 
of the critical temperature and the scalar operator on the presence of the dark matter sector. 

To commence with, one considers the background of five-dimensional black hole given by the line element 

J.2 

ds^ = -g{r) dt^ H— -p-- + —r {dx^ + dy^ + dz'^), (39) 

g{r) 


where g{r) = r'^jLp' — The Hawking temperature for the black hole has the form Tbh 

z-coordinate the equations of motion imply 


V'" 


T] = 




q2 02 


9 zj \ g 

1 2 0 ^ 0 rl 

z a q 

di z a , 
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= 0 , 


0 = 0, 


r+/7r. In 


(40) 

(41) 

(42) 


where the prime denotes the derivation with respect to z-coordinate. In the following, as in the preceding sections, 
we set L = 1. In order to solve the above equations we need to impose the adequate boundary conditions. At the 
black hole horizon z = 1, it is required that 0(1) = 0 and 0(1) should be finite. The first requirement is needed for 
the [/(l)-gauge field to have the finite form, the second one exhibits that the black hole has a scalar hair on the event 
horizon. 

When the temperature T tends to the critical value Tc from below, the condensation approaches zero, 0—^0. In 
this limit we write the equation for 0 field as 


1 

z 


0 ' ~ 0 . 


(43) 
















Its general solution is of the form 4>{z) = ci + C 2 Z^, which together with the aforementioned boundary condition at 
the horizon leads to (j){z) = Ci(l — z^) . Then, having in mind that near the boundary of the bulk the fields behave as 




^ pz'^/rl, 

^A_ J.A+ ’ 


(44) 

(45) 


one arrives at the conclusion that near the critical temperature the At gauge field component will behave as ~ 
A r+ (1 — z^), where we have denoted by A = p/r^. p and p have the same interpretation as in the preceding sections, 
i.e., denote the chemical potential and the charge density, respectively. Moreover, A± have the same values as in 
Sec.II, but in order to distinguish the nature of the phase transition, we set the new notation for the old quantity. 
In what follows we concentrate on the A = A+ and the masses fulfilling the Breitenlohner-Freedman spacetime 
stability condition. It is important to note that we are looking for the parameter A which defines the value of both p 
and p. It depends on the temperature given by T = Tbh = 

On the other hand, the scalar field '0 can be cast near the boundary in the form 

0 U^o- < C > ^ G{z), (46) 


where, we set G(0) = 1, G'(0) = 0. Inserting this expression into the equation (HOI) and using 0 = A r_|_ (1 — z^) one 
finds the scalar held equation of motion, which can be easily rewritten as the Sturm-Liouvile eigenvalue problem 

{Viz) G'iz))' - Qiz) G{z) + A^ n{z) G{z) = 0, (47) 


where one has dehned the quantities 

V = 5, (48) 

Q = -A(A -l)g z^^-^ - A g + m^rl z^^-\ (49) 

^2A-5 

= - ri (1 - zY q^. (50) 
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The minimum value of A^ can be found from variation of the functional given by 

^2 _ Jo [G'jz? 'Pjz) + Q{z) G^(z)] 
lodz n{z) G^iz) 

where the trial function is assumed to be given by G(z) = 1 — a z^. Note that in the variable z, the function 
g{z) = r^(l/z^ —z^) and the value of event horizon radius r_|_ factors out of the expression for A^ making it independent 
on Tbh- Because of that and the fact that the analysis is valid close to the transition temperature Tc ^ Tbh, one 
readily hnds that 


T, = 


7r3 A„ 


1/3 


(52) 


It is important to note that the value Xmin results from the variation of Sturm-Liouville functional (15111 . Holographic 
superconductors with different values of p are characterized by different transition temperatures. This reminds the 
charge carrier concentration dependence of the local pairing superconductors [4l|. However, the charge density 
dependence of the critical temperature Tc oc p^^^ found here markedly differs from the known dependence Tc oc 
for the Bose-Einstein condensation of the low density n superconductors with local pairs of charged hard core bosons 
[4l[. Whether the difference is the hallmark of the strong coupling behavior remains to be seen. In fact it has been 
suggested @ that the condensation transition in holographic models is closer to the Bose-Einstein condensation than 
BCS-like symmetry breaking phase transition. It has to be reminded that the critical temperature does not depend 
on the carrier concentration n in the standard weak c oup ling superconductors described by the BCS theory. Instead 
it depends on the density of states at the Fermi level . 

It can be observed that the critical temperature does not depend on the dark matter sector. This conclusion is in 
accord with our previous studies (sH , where it was revealed that the backreaction effects introduce the dependence of 
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Tc on the dark matter sector. In the case under consideration we restrict our investigations to the probe limit case, 
therefore no influence is spotted. 

In the Table IIIII we have presented the results of the calculation of the superconducting transition temperature 
assuming the charge density p = 1. The presented theory is valid for arbitrary allowed values of m^, but for 
illustration we have chosen some exemplary values of it and calculated and the critical temperature T^. Both 
of these parameters are presented in the Table IIIII together with the value of the parameter a, which minimizes the 
functional (IFTl) . 


A. Condensation values 

In this subsection our main task will be to find the influence of the dark matter sector field on the condensation 
operator. Near the critical temperature, the equation for gauge field can be rewritten as 


2 r? 


,2A-4 


a g 


G\z) A, 


(53) 


where we have denoted A =< C Having in mind the fact that near the critical temperature the A quantity 

is small, one expands cj) near z —>• 0. On this account, we can write the following: 


— = A (1-^2)+AIx(z) + .... 
r+ 

In the next step, comparing the coefficients in z^-order terms, we reveal that 

On the other hand, considering the relation (IMl) and the equation of motion for cj) field, we get 

// x' “2, \ 2n „2A-4 ,-i2/ 


^ z 


<5 / 


(1-z") G\z). 


Consequently, proceeding as in the last section, it can be verified that the following is satisfied 


X (0) —- U-)-o— —2 A 


« / 


(54) 


(55) 


(56) 


(57) 


On evaluating the expression for A, in the case when T —> the condensation operator in question is provided by 


<c>= 


B 


TcY 


l-^= 

n 


'a <C> 


no dark sect: 


(58) 


where by < C >no dark sect we denoted the value in the theory without dark matter sector. The term B is given by 
the relation 


B = 2 


A 2 

Jo 9 


(59) 


One can see that the condensation operator depends on the a constant coupling of the dark matter sector. The 
bigger is the a-coupling, the easier condensation forms. From the point of view of the AdS/CFT correspondence, 
the operator in question can be interpreted as the operator for pairing mechanism. The bigger expectation value it 
achieves, the harder condensation occurs. In order to better understand the dependence of the condensate on the 
dark matter sector we factor out its dependence on the rest of the parameters in the following way: 


^ <C> _ /rzL T 

< O >norm— = V OiG y ^ ~ 


(60) 


where C = \ and < G >r 


represent renormalized value of the condensate. The factor C depends on all 

the remaining parameters (except temperature) and its typical values are presented In Table IIIII These facts can 
potentially constitute the way of determining the dark matter sector in future ’possible’ superconductor experiments. 

We remark that, in the case of Gauss-Bonnet theory [i^, the value of < C > Is dependent on the higher curvature 
term corrections. When It grows, the value of the operator also increases. This conclusion is in agreement with the 
previous studies (see, e.g., (39| and references therein). 
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A 


A^ 

dmin 

1 — 1 

II 

Q. 

C 

5 

2 

-15 

4 

9.586 

0.619 

0.218 

55.67 

6 

2 

12 

4 

18,22 

0.721 

0.196 

137.8 

7 

2 

7 

4 

30.50 

0.797 

0.180 

326.7 

8 

_2_ 

0 

46.89 

0.853 

0.168 

748.3 


TABLE III: Calculated values of the prefactors of minimizing the functional (ISTI) . the superconducting transition 
temperature of the superconductor with charge denisity p = 1 and the prefactor C defined in Eq. (1601) . In all the 
above examples we assumed g = I.O and selected the trial function as G = 1 — az^. We also provide the parameter a 

minimizing the functional (1511) . 


IV. HOLOGRAPHIC DROPLET IN S-WAVE INSULATOR/SUPERCONDUCTOR PHASE 

TRANSITION 

The term superconducting droplet refers to the solutions that are confined in space and rapidly decay at large 
distances. This happens if the studied superconductor is exposed to strong external magnetic field. The size of the 
conhning region diminishes with the increase of the magnetic field. In this section we shall investigate the onset of 
the transition by studying the marginally stable modes of scalar perturbations. It has been shown earlier [l^, 113 

that in the holographic approach the marginally stable modes signal the appearance of the insulator-superconductor 
transition. 

Here we are interested in the insulator - superconductor phase transition in the AdS solitonic background coupled 
to the dark matter sector. The presence of the magnetic field introduced via dark matter potential enables us to 
examine the droplet solution via the aforementioned technique. 

It turned out that the quasi-normal modes (QNMs) technique occurred as a method of examining stability of a 
spacetime background (l^ . In the case when the imaginary part of QNMs is negative, the modes decrease in time 
and result in disappearance of perturbations (background is stable against perturbations). On the other hand, when 
imaginary part is positive, the background is unstable against the perturbations in question. The marginally stable 
modes are the modes which frequencies go to zero (w = 0) near the critical point of the phase transition. Their 
appearance thus signals the phase transitions [HI, . 

To commence with, we consider AdS soliton metric, taking into account symmetry of the problem in question and 
rewrite the line element used in Sec.II in the following coordinates (t, , r, p, v?, 0) as 

dt^ -I- -I- /(r) dtp'^ + (dp^ -|- dd^), (61) 

f(r) 

Eurther, we assume the existence, in addition to a constant chemical potential bounded with At gauge field component, 
the Bg potential which corresponds to the dark matter gauge field and is proportional to the constant value of magnetic 
field B 


At = p, Bg = ^B p^. (62) 

The above ansatz stems from the fact that one considers gauge sector close to the critical point of the phase transition, 
i.e., p Pc and '0 ^ 0, as well as, the implementation of the polar coordinates in order to envisage the symmetry of 
the problem in question [Hj |. 

Having in mind the exact form of Bg, given by the above relation, we calculate from equation of motion 

+ I = 0, (63) 


the Ag component, which yields 


Ag = 


Di r p^ 


B 




(64) 


Di and D2 are integration constants. To proceed further, we assume Ag to be the function of p only, which implies 
that £>1 and D 2 have to be equal zero. Our studies are devoted to the probe limit, i.e., the C/(I)-gauge fields and 
scalar one do not backreact on the AdS soliton background metric. Without the condensate (i.e. for ■0 = 0) the 
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solution of the equations of motion for the A gauge field components are given by the Ag = and At = p. We 

are interested in finding the solution of if) equation close to the critical chemical potential ^ ~ /ic, where the value of 
the scalar field reaches nearly zero. It implies that one can treat ip field as a probe into the background consisting of 
AdS Schwarzschild soliton with constant electric and magnetic field. 

The equation describing ip field yields 


Vfj,y^ip — Afj^A^ Ip — ip = 0. 


By virtue of the above, the explicit form of the equation for ip field may be written as 


d^ip 


1 2 , 

— + -jdrip + j^d^ip- 


1 


1 


^2 f 


J.2 f 

2 2 2 2 

q p — m r - - - 


d'^ip 


1 1 


f p 


dp[p dpip) 


Ip = 0. 


In order to solve the above equation we choose an ansatz for ip field 

Ip = F{r, t) H{(p) U{p). 


(65) 


( 66 ) 


(67) 


This form enables us to separate variables. After a simple algebra we arrive at the following set of equations 


dlF + 


3 

r f 


drF - d?F + ^ 


p2 f 


p2 f 


2 2 2 2 
q — m r — 


-k^] F = 0, 


~ir - 


dp p d-pU - 


q^ B'^ p^ 


U = -k^ U, 


( 68 ) 

(69) 

(70) 


From periodicity property H{ip) = F[{ip tt FjrQ) of H{ip) we identify that A = 2 ro n/L, where n G Z. In what 
follows without loss of the generality we set tq = 1 and L = 1 what leads to A = 2n. We expect that the lowest mode 
will be first to condense and result in the most stable solution. 

The equation for U{p) is a two-dimensional harmonic oscillator one. In order to solve it we recall that the function 
U{p) should satisfy the boundary conditions C/(p —>■ oo) = 0. It is possible to investigate such kind of differential 
equation by Frobenius method. The customary procedure is first to factor out the behaviour of the relevant solutions 
at infinity by setting 


Uip) = D{p\ (71) 

which results in Hermite’a type of equation. Inserting dzlD into the underlying equation, we obtain the expected type 
of the differential solution if the condition 


A" = 



(72) 


is satisfied. Then, the resulting equation yields 


\ dpi p dpD ) - 2 A p dpD + {k^ -2 A)D = 0. 


(73) 


In order to find the exact form of D{p), one sets D{p) = EB)- 

Note that D{p) = const is the well know lowest energy solution of the harmonic oscillator. It leads to the condition 


= 2 A =1 I . (74) 

a 

Equation GD) shows that in the presence of ma gne tic field the superconducting region is confined in space. For the 
chosen solution D{p) = const it forms a droplet of radius 


J dp pU if) ^ 1 ^ / a 

J dp U (p) v^27riV y 4: TT q B 


( 75 ) 
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In the field theory it corresponds to the condensation in the lowest Landau level . As far as the time dependence 
of the F(r, t) is concerned, we substitute it in the form F(r, t) = ^(^)- The requirement concerning marginally 

stable modes leads to the condition w = 0. Redefining the coordinates as z = r+/r enables to arrive at the equation 
given by 


9zR{z) + ( — - j dzR{z) + ^2 


2 2 2 R g m2 4 n^\ 

,2,y‘l l‘ - —-TT-jjyJ «W = »- 


To solve it close to when the field ijj fa 0 we introduce a correction function 0(z) in the form 

R{z) <0,> z^* 0(z), 


(76) 


(77) 


with the boundary conditions 0(0) = 1 and 0'(O) = 0. After some algebra, the resulting equation can be converted 
into the standard Sturm-Liouville eigenvalue equation, which can be rewritten as 


dz 



b{z) e + S^ ciz) 0 = 0, 


(78) 


where — 2 q B/a and the remaining quantities are defined by the relations 

a{z) = f 


b{z) = -f 




/ 


z^f 


4 


m 


_2Ai-3 


c{z) = z- 

The eigenvalues of 6'^ can be found by the method of minimizing the functional 

2Bq f,^dz(e'(zfa(z) + b(z)e(z)^) 


ir2 11 

(5 =q 


Q 


fo dz (c(z) 02(z)) 


(79) 

(80) 
(81) 


(82) 


In order to estimate 5^, we choose function 0 ( 2 ) = 1 — a z^. Minimization of the functional provides an estimation 
of the value of S which depends on and n, resulting from the periodicity property of F{{ip). The above relation 
between S^, critical chemical potential can be rewritten as 


Me = 



(83) 


It follows that Me depends on the coupling to the dark matter sector a. Interestingly, for constant magnetic field B 
the critical chemical potential diverges for a —>■ 0. It means that at the constant magnetic field the condensation is 
harder to occur for smaller values of a. However, one has to remember that the zero value of a is not allowed, as 
simultaneously one has to take i? = 0, so without dark matter field the standard relation, valid for zero magnetic field 
[ 2 ^ <52 = is recovered. The increase of the magnetic field causes the increase of Mc, which in turn eventuates 
in the harder condensation. The aforementioned behaviour is depicted in Figl2]for two values of the magnetic field. 
The discussed increase of Mc for small a is clearly visible both for H = 0.1 and for B = 1. 


V. DISCUSSION AND CONCLUSION 

The main aim of our paper is to find the quantitative or at least qualitative imprints of the dark matter sector 
on the properties of s-wave holographic superconductor phase transitions. The unordinary features might constitute 
the possible hints for future experiments testing the considered model of dark matter. In the model in question, 
apart from the electromagnetic matter field we have taken into account the dark matter sector described by another 
[/(l)-gauge field, bounded with the Maxwell field by the coupling constant a. 

The models where dark matter is a part of a larger sector which interacts with visible matter were successesfully 
implemented as the possible explanations of various astrophysical anomalous observations like the excess of electrons 
in Galaxy having energies of a few GeV and TeV, gamma rays of 511 keV [55|. There were also efforts to find new 

a sics explaining the anomalous muon magnetic moment, possible implication for parity violation, rare meson decays 
, as well as, to provide some implications of boson and dark boson mixing for high energy experiments This 
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problem is of a great importance especially in the light of the latest claim of nongravitational interactions of dark 
matter in colliding galaxy clusters [^, which can disfavor some extensions of the Standard Model. 

In the paper we have discussed analytically various phase transitions toward the s-wave holographic superconductor 
in the probe limit. Our results for a = 0 agree with the previous numerical and analytical studies of holographic 

superconductor transitions. The coupling between ordinary and dark matter changes the values of the parameters at 
the transition, however, only quantitatively for s-wave superconductors. 

To make contact between transitions studied here and those known in the condensed matter systems, let us re¬ 
call some basic facts from the latter field of research. In the condensed matter systems there exist a number of 
metal/insulator transitions. They differ by the role played by the interactions between carriers and the lack of period¬ 
icity of the underlying crystal lattice. The metal/insulator transition may appear when the carriers strongly interact 
with each other. On the physical grounds it can be argued that in such strongly correlated system the electron 
movement is hindered by the repulsive interactions and as a result, the insulator may form. The literature on the 
Mott-Hubbard metal/insulator transition [lailil is vast and the transition itself is still not fully understood [s^. 
The other interesting transition appears in the system which is not translationally invariant. The transition is driven 
by disorder. Strong disorder makes some states localized and thus unable to carry the electric current. This is called 
Anderson metal/insulator transition 15311. In real materials one usually finds transitions in which both interactions 
and disorder do play an essential role |54| . 

Our study is restricted to the system which is periodic on the boundary and the insulator/metal transition we 
study should be related to the Mott-Hubbard one. The sequel of the transitions we are discussing, i.e., insula¬ 
tor/superconductor at zero temperature and metal/superconductor at higher temperatures are realized in high temper¬ 
ature cuprate and iron superconductors with the increase of the carrier doping. The T = 0 insulator/superconductor 
transition is an analog of the Hawking - Page like soliton-black hole transition. In real superconductors increase of 
charge density beyond the upper limit induces (at low temperatures) a reverse superconductor to metal transition, 
which seemingly has not been hitherto found in holographic analogy. 

Treating as a spectral parameter we analyze the behaviour of it in s-wave insulator/superconductor phase 
transition. We did not observe the influence of the dark matter sector on this quantity. The charge density is 
proportional to the difference — ^c), and also does not depend on the dark matter sector in the probe limit. In the 
case of the scalar operator it was revealed that it is proportional to and represents the second order phase 

transition. The critical exponent of the considered system has the mean-held value, while the proportionality factor 
is subject to dark matter coupling constant dependence. The smaller value of a is considered, the greater factor one 
obtains, i.e., the harder condensation happens (it will also be the case in s-wave holographic metal/superconductor 
phase transition). 

The same conclusions were drawn studying s-wave type of the transition in question in Gauss-Bonnet theory. The 
Gauss-Bonnet coupling which envisages the inhuence of higher curvature corrections, does inhuence the considered 
factor and the form of linear dependence survives. 

In the case of holographic metal/superconductor phase transition, the critical temperature does not depend on dark 
matter sector. The conclusion is in accord with our previous studies [3^ . where it was shown that the backreaction 
effects introduced the dependence of Tc on dark matter sector. On the contrary, the condensation operator reveals 
a linear dependence on the a-coupling constant. The bigger dark matter coupling constant one considers, the easier 
condensation forms. 

Examining s-wave droplet insulator/superconductor phase transition, it was found that the chemical potential and 
magnetic field were bounded with the linear dependence. Not only does the magnetic held inhuence the condensation 
but also the coupling constant of dark matter sector does this. The increase of the magnetic held causes the increase 
of /ic, which eventuates in the harder condensation. For the hxed value of the magnetic held, it happens that the 
smaller value of the dark matter sector coupling constant one chooses, the harder condensation takes place. A very 
similar behaviour was envisaged in the holographic droplet in p-wave insulator/superconductor phase transition case. 

To conclude, we remark that there are some points which are in contrast to the ordinary behaviour (without dark 
matter sector) during the aforementioned phase transitions, which may constitute indicators for the future experiments 
for detecting dark matter and elucidating its nature. Testing s-wave holographic superconductors with dark matter 
sector is the only tip of the iceberg and some more complicated models like p-wave or -f ipy should be taken 
into account. We hope to investigate these problems elsewhere. Our preliminary calculations show the stronger 
modifications of the p-wave superconductor characteristics by the dark matter coupling a. 


Appendix A: Integration constants in Eq. (1321) 


Equation ()32ll provides formal solution of the differential equation (1^^ and contains two integration constants. The 
value of ^(0) enters the prefactors T and B of the dependence of the condensation operator < Oi> and the density 
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p on the chemical potential. 

To find numerical values of the coefficients T and B we note that the integrals entering equation (15^ still contain 
terms singular at z = 1 which have to be eliminated by the proper choice of constants. To this end we take 


B(y) = 1 — ay^ and for numerical evaluation of the constants use the values of a = 
functional (l20ll for p,^. Evaluating the (non-singular) integral over y in Ea. (l32l) leads to 

1 


which minimize the 


9^ rl dy y 


2\i-3 7^2 


F\y)=q^ 


2Ai 




A, 2(A, - l)y2 2(A, + 1) 


2 2 
= 9 ’’o 


.2Ai 


a 


1 


Ai 2{\j - l)z2 2(A, + 1)) _ 


Xi 2(Aj — 1) 2(Ai + 1)/ 


Inserting the above result denoted as R{z) into (l32l) and performing the integral over z one gets 


(Al) 


^(0) = Cl - 


dz- 


[C 2 + R{z)] = Cl - [^^(log(l -k Z^) - log(l - z^)) + rlw{z) 


1 - 


J 0 


where we have denoted 
Wiz) = 


2Ai(l -l- a^(Ai — 1) -l- Ai)z^^’~'’‘^ 2A'i(l, ^ ^ — (2 -I- Ai)(—2z^^* (A^ -I- 1) 


-f4a(A, - 2-F'i(1, + (a^A,(Ai - 1) - Ai(A 7 -k 1) - 2a(A- - 1)) log(l - z^) 

-|-(— 2a -l- Ai -l- a^ Xi -I- A^ -I- 2aX^ — a^A^) log(l -I- 2 ^)) /[8(Ai — l)Ai(Ai -I- l)(Ai -I- 2)]. 


(A2) 


(A3) 


Here 2 ^ 1 ( 0 , b; c, z) is the hypergeometric function |4dl| . For the special values of parameters with c = a + 6, as in the 
above expression, it diverges for z —>■ 1 and takes the form [i^ 


lini 2 -Fi(a, b;a + b,z) = - log(l - z). 

i (a)i (0) 

It is easy to check that IT(0) = 0. On the other hand the boundary condition ^(1) = 0 requires 


C2 = 


2 2 

9 eg 

A-1’ 


and 


2 2 

Cl = -9 Cq 


log 2 


_4(Ai — 1) 4Ai(Ai — 1) \4Ai(Ai + 2) 4(Ai + l)(Ai + 2) 8(Ai — l)(Ai + 2) 
Introducing the obtained results into the eouation d^Sl) . one gets the required formula 

< Oi >= T^ 


log 2 


(A4) 


(A5) 


(A6) 


(A7) 


where the prefactor T = (a-i)^ log 2 ■ other hand, the charge density yields 

4(Ai - I)/o^y^^’“^F2(i/)dy 
-(a-l)Uog2- 


(AS) 


Another way to find the required value ^(0) is by direct solution of the equation (l29l) . For general value of the 
parameter A^ it is given in terms of the hypergeometric 27 ^ 1 ( 0 , 6; c, z) functions as 


= C 2 + 


2 2 
9 


4(Ai + 2)(Ai + I)Ai(Ai — I) 


A.(a2(A. - I) + Al + 1)z^^^+SFi{1, ^ + I; ^ + 2, z^) 


— {Xi + 2) ( 2a{Xi — I)z^^*~*'^ 2 Ti(I, —-—; —-—, z^) + {Xi + I)(—+ Ai(Ai — I)Ci log 


I - z" 

TT^ 


(A9) 


Remembering that [i^ lim^-n 2 -Fi(l, ^ + F, ^ + 2, z^) = —(1 + ^) log(I — z^) and lim^^o 2 A'i(I, ^ +1; ^ -I- 2, z'^) = I 
one chooses the constants Ci and C 2 in such a way that the solution is finite with x(I) = 0 and finds ^(0) = C 2 with 
C 2 given by the formula (IA6I) above. 
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FIG. 1: (color online). The dependence of the prefactor of the condensate function on the coupling constant a. The 
trial function was chosen as F = 1 — az^ and the charge was set equal to g = 1. 




FIG. 2: (color online). The dependence of the critical chemical potential on a in the droplet case, for the fixed value 
of the magnetic field. We set B = 0.1 (for the left panel) and B = 1 (for the right panel), the rest of the parameters 
are equal to g = 1.0, n = 0, = —^(A = |), = —3 (A = 3). 














